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Abstract 

We compute exact solutions of two-matrix models, i.e. detailed genus by genus expressions 
for the correlation functions of these theories, calculated without any approximation. We distin- 
guish between two types of models, the unconstrained and the constrained ones. Unconstrained 
two-matrix models represent perturbations of c = 1 string theory, while the constrained ones 
correspond to topological field theories coupled to topological gravity. Among the latter we 
treat in particular detail the ones based on the KdV and on the Boussinesq hierarchies. 
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1 Introduction 



Matrix models represent sums over discretizations of Riemann surfaces, possibly with some 
additional interactions. They are believed to provide a (discrete) description of two dimensional 
gravity coupled to matter. One-matrix models have been widely investigated, but their content 
is rather poor. The structure of multi-matrix models is much richer but not yet known as 
carefully as for one-matrix models (see |], |],[|, §,§,§,§, i,@,0,O,[0],0,0). 

In this paper we concentrate on two-matrix models with bilinear coupling and show how to 
find exact solutions. By exact solutions we mean detailed genus by genus expressions for the 
correlation functions of these theories, calculated without any approximation. In particular, we 
do not limit ourselves to exhibiting recursion relations which allow one to compute correlators, 
but develop techniques to explicitly solve them. 



The idea at the basis of our treatment of two-matrix models, outlined in |15|, is to transform 
the initial functional integral problem into a discrete integrable linear system subject to some 
constraints (the coupling conditions). We end up in this way with the discrete Toda lattice 
hierarchy. The latter underlies all our calculations: our aim is to compute the correlation 
functions (CF's) of each model, which in turn may be expressed in terms of the integrable 
flows of the Toda hierarchy, subject to the coupling conditions. This is the general setting for 
unconstrained two-matrix models (which simply denote two-matrix models defined by specific 
potentials without any further conditions). 

We will also consider other models, obtained by suitably constraining the previous (uncon- 
strained) models. Their correlators can be expressed either in terms of the flows of a reduced 
differential hierarchy or in terms of suitably renormalized flows of the discrete Toda hierarchy. 
In the process of solving these models we find a new way of extracting integrable differential 
hierarchies from the Toda lattice hierarchies. 

Unconstrained two-matrix models describe various perturbations of c = 1 string theory at 
the self-dual point. Constrained models correspond to well-known topological field theories 
coupled to topological gravity. 

The paper is organized as follows. In section 2 we collect the results obtained in previous 
papers which will be necessary in the following. In section 3 we discuss and calculate CF's of 
unconstrained models. Section 4 is a short summary of how to obtain differential hierarchies 
from the Toda lattice flows and reduced hierarchies via hamiltonian reduction. In section 5, 
from the reduced hierarchies we construct a series of models, named reduced models, which have 
a topological field theory meaning: we show in particular how to compute all genus correlators. 
The reduced models turn out to be imbedded into the constrained two-matrix models, which 
are studied and solved in section 6. 

2 General properties of two— matrix models. 

The model of two Hermitean N x N matrices Mi and M2, is introduced in terms of the partition 
function 

Zjv(i, c) = J dMidMae*''^, U = Vi + V2 + gM^Mi (2.1) 

with potentials 

Pa 

^a = EW^a a = 1,2. (2.2) 

r=l 
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where pa are finite numbers. These potentials define the model. We denote by A^pi,p2 the 
corresponding two-matrix model. 

We are interested in computing correlation functions (CF's) of the operators 

Tfc = trMf , cjfc = trM^, V/c, x = tr{MiM2) 

For this reason we complete the above model by replacing (^]^) with the more general potentials 



^a = ^W^a> a = 1,2 (2.3) 



r=l 

where ta,r = ta,r for r < Pa- The CF's are defined by 

<Tr,... Tr„as, ■ ■ ■ CT^^x' >= JT. TT. TT. TT. TnlnZN{t, q) (2.4) 

where, in the RHS, all the ta,r except ta,r are set to zero. 

In other words we have embedded the original couplings ta,r into two infinite sets of couplings. 
Therefore we have two types of couplings. The first type consists of those couplings (the barred 
ones) that define the model: they represent the true dynamical parameters of the theory; they 
are kept non-vanishing throughout the calculations. The second type encompasses the remaining 
couplings, which are introduced only for computational purposes and are set to zero in formulas 



like (2^). In terms of ordinary field theory the former are analogous to the interaction couplings, 
while the latter correspond to external sources (coupled to composite operators). From now on 
we will not make any formal distinction between them. Case by case we will specify which are 
the interaction couplings and which are the external ones. Finally, it is sometime convenient to 
consider on the same footing as the couplings and to set ii^o = ^2,o = 

The path integral (|2.1| ) is an ordinary integral in the matrix entries and it is certainly well 
defined as long as a negative coupling fa,r with highest even r guarantees that the measure is 
square-integrable and decreases more than polynomially at infinity. For the time being let us 
suppose that this is so. Later on we will extend our problem to a larger coupling space. 

We briefly recall the ordinary procedure to calculate the partition function. It consists of 
three steps [17|,[18|,[|l^: (i) one integrates out the angular part so that only the integrations 



over the eigenvalues are left; (ii) one introduces the orthogonal monic polynomials 

(,n{^i) = A" + lower powers, '>]ni^2) = A2 + lower powers 

which satisfy the orthogonality relations 

J c?AidA2en(Ai)e^(^i'^2)r?„(A2) = hn{t, c)5nm (2.5) 

where 

^(Ai, A2) = yi(Ai) + F2(A2) + CA1A2 



(iii), using the orthogonality relation ( p.5[ ) and the properties of the Vandermonde determi- 
nants, one can easily calculate the partition function 

ZNit,c) = const m Y[ hi (2.6) 

i=0 



3 



2.1 From path integral to integrable systems 



From ( ^.61 ) we see that knowing the partition function means knowing the coefficients hn{t, c). 

The crucial point, from our point of view, is that the information concerning the latter can 
be encoded in 1) a suitable linear system subject to certain 2) coupling conditions, together 
with 3) relations that allows us to reconstruct Z^. But before we pass to these three elements 
we need some convenient notations. For any matrix M, we define 

M = H-'MH, H,j=hi6ij, Mij = Mji, Mi{j) = Mj^^.i. 

As usual we introduce the natural gradation 

deg[Eij] = j - i, where {Ei,j)k,i = ^i,k^j,i 

and, for any given matrix M, if all its non-zero elements have degrees in the interval [a, 6], then 
we will simply write: M G [a, Moreover M+ will denote the upper triangular part of M 
(including the main diagonal), while M_ = M — M_|_. We will write 

N-l 

TV(M) = Mu 

1=0 

The latter operation will be referred to as taking the finite trace. 

Next we pass from the basis of orthogonal polynomials to the basis of orthogonal functions 

^n(Ai) = e^i(^i)en(Ai), <I>„(A2) = e^^(^^)r/„(A2). 

The orthogonality relation (^]^) becomes 

d\id\2^n{\l)e'^'^^^ra{X2) = SnmKit, c). (2.7) 



We will denote by ^ the semi-infinite column vector (^fg, ^'i, ^'2) •••))* and by $ the vector 
($0, ^1,^2, •••,)*• 

Then we introduce the following Q-type matrices 

j d\id\2^n{\l)Ke^^'^^^m{\2) = Qnm{a)h^ = Qmn{a)K, a = 1,2. (2.8) 

Both Q{1) and Q(2) are Jacobi matrices: their pure upper triangular part is /+ = J2i^i,i+i- 
Beside the above Q matrices, we will need two P-type matrices, defined by 

dAi(iA2(^^'n(Ai))e'=^^^2$™(A2) = Pnm{l)hm (2.9) 

d 

For later use we also introduce 

d 



j dXidX2^n{Xl)e'^'^'{-Qj^'^m{X2))=Pn,n{2)hn (2.10) 



J dAidA2(^en(Ai))e^^(^^)+^^(^^)+^^^^^r?„(A2) = P°^(l)/i^ (2.11) 
dAidA2Cn(Ai)e^i(^^)+^^(^^)+^^^^^(^r?„(A2))=P°„,(2)/i„ (2.12) 
Let us come now to the three elements announced above. 
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1) Coupling conditions. The two matrices ( |2.8D we introduced above are not completely 
independent. More precisely both (5(a) 's can be expressed in terms of only one of them and one 
matrix P. Expressing the trivial fact that the integral of the total derivative of the integrand 
in eq.(2.7) with respect to Ai and A2 vanishes, we can easily derive the constraints or coupling 
conditions 

P(l)+cg(2) =0, cQ(l)+P(2) = 0, (2.13) 

Prom the coupling conditions it follows at once that, if we set to zero the external couplings, 

Q(q) G [-m„,n„], a = 1,2 



where 



and 



mi=p2- 1, 



ni = 1, 



7712 



n2 = Pi 



1 



where pa, ct = 1,2 is the highest order of the interacting part of the potential (see above). 

2) The associated linear systems. The derivation of the linear systems associated to our 
matrix model is very simple. We take the derivatives of eqs.(2.7) with respect to the time 
parameters ta,T, and use eqs.( |2.8| ). We get in this way the time evolution of or discrete linear 
system I: 

^ ^(Ai) = Q^(l)+^(Ai), 
-Q'^(2)_vi;(Ai), 



a 



(2.14) 



5*2 



^(Ai) 



|^vI;(Ai)=P(l)^(Ai). 



The corresponding consistency conditions are 

1 



-^g(i) = [Q(i), qH») 

Ota.k 



a 



1,2 



(2.15a) 
(2.15b) 



In a similar way we can get the time evolution of ^ via a discrete linear system II, whose 
consistency conditions are 



[2(2), P(2)] = l, 



d 



a,k 



-Q(2) = [Q'=(a)+, Q(2)] 



(2.16a) 
(2.16b) 



We recall that one can write down flows for P{1) and P{2), but they will not be used in this 
paper. 

3) Reconstruction f ormulae. The third element announced above is the link between the 
quantities that appear in the linear system and in the coupling conditions with the original 
partition function. We have 



d 



do 



■lnZAr(t,c) = Tr(Q'^(a)), a = 1,2 



(2.17) 
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It is evident that, by using the equations ( 2.15b , ^16b ) above we can express all the derivatives 
of Zi\f in terms of the elements of the Q matrices. For example 



dtiAdto 



-lnZ7v(t,c) = (q^'ia 



N,N-l 



a 



1,2 



(2.18) 



and so on. We recall that the derivatives of F{N,t, c) = \nZN{t, c) are nothing but the correla- 
tion functions of the model. 

We can summarize the content of this section in the following 

Proposition 2.1. The correlators (2.4) can be expressed in terms of the entries of the 



matrices Q{1) and Q{2) via eg. ( 2. 11 ) and the like. In turn, these matrices must satisfy the 
coupling conditions ( ^.13^ and the consistency conditions ( 2. 15q - 2. 16b ). 

Knowing all the derivatives with respect to the coupling parameters we can reconstruct the 
partition function up to an overall integration constant (depending only on A^). The recon- 
structed free energy F will be a power series in the external couplings. 



This theorem was proven in |15|. It is a rigorous result when, for example, highest negative 
even couplings guarantee that the measure in ( |2.5| ) is square-integrable and decreases more 
then polynomially at infinity. But for generic values of the couplings the above derivations are 
merely heuristic. 

However we notice that the consistency and coupling conditions make sense for any value of 
the couplings, and also when the couplings are infinite in number. In the latter case eqs.( p.l5b ) 
and ( 2.16b| ) form nothing but a very well-known discrete integrable hierarchy, the Toda lattice 
hierarchy, ||2^ (see also |^l|). Prom these considerations it is clearly very convenient to refer to 
the integrable system formulation rather then to the original path integral formulation of our 
problem. This allows us not only to extend our problem to a larger region of the parameter 
space, but also to make full use of integrability. Therefore we shift from the original problem 
to the new formulation: 

We call (unconstrained) two-matrix models all the models obtained by specifying a partition 
of the couplings between internal and external couplings. Each such model is based on the Toda 
lattice hierarchy and characterized by specific coupling conditions. 

Statement of the problem. Solve the integrable Toda lattice hierarchy subject to the 
coupling conditions specific of a particular unconstrained model and compute the correlators as 
functions of the internal couplings via the relations 



< Tr >= Tr Q(l 



<as >=Tr Q(2) 



(2.19) 



and the like. 

Once all the correlators are known, one can reconstruct f\ the free energy F by means of 



d 



dti 



■F{N,t,c) =< Tr >, 



d 



dt 



2,r 



■F{N,t,c) =< ar > 



(2.20) 



F will be a formal power series in the infinite set of external couplings. 
Henceforth this will be the setup we refer to. 

To end this subsection, we collect a few formulas we will need later on. First, we will be 
using the following coordinatization of the Jacobi matrices 



mi 



"12 



Q(l) =/+ + ^^a,(i)i^i,i_,, 

i 1=0 



Q{2) = 1+ + Y,Y.^i{i)Ei^i_i (2.21) 

i 1=0 



*Up to a constant depending only on TV. There is a way to determine this constant too, see [Q, but we will not 
dwell upon this point here. 
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One can immediately see that 



(q+(1)) ..= 6j,i+i + aoii)5^,J, (q_{2))^= (2.22) 



where R{i + 1) = hi^i/hi. As a consequence of this coordinatization, eq.( p.l8 ) gives in partic- 
ular the two important relations 

^F{N,t,c)=ai{N), (2.23) 
"^^1,1 

Finally we write down explicitly the and t2, inflows, which will play a very important role 
in what follows 

^^oKi) = a/+i(j + 1) - ai+iii) + a/(j)(«o(j) - ao(j - O) (2.24a) 
^ ai{j) = R{j - I + l)ai^i{j) - R{j)ai^i{j - 1) (2.24b) 



5^2,1 

d 



dt 



2.2 Woo constraints. 



kij) = R{j - I + l)6,_i(j) - Rij)bi_,{j - 1) (2.24c) 
^ -6z(j) = bi+iU + 1) - 6;+i(j) + biU){boij) - boU - I)) (2.24d) 



To solve the above stated problem we have to solve the flow equations of the Toda lattice 
hierarchy subject to the coupling conditions (|2.13D . There is a way to put together flow equations 
and coupling conditions that lead to an elegant algebraic structure, the W constraints: 

Proposition 2.2 The partition function of the unconstrained two-matrix models satisfies 
the following W constraints: 

Wj{^ZN{t,c)=0, Wj{^ZN{t,c)=0 r>0; n>-r, (2.25) 

where 

lyM ^ (-c)"£M(i)_£b;+"l(2) (2.26a) 

lyM ^ (-c)"£M(2) -£[!+"'(!) (2.26b) 

The generators £n^(l) are differential operators involving and ti^fc, while £n'(2) have 



the same form with ti ^ replaced by t2^k- One of the remarkable aspects of (2.25) is that the 



dependence on the coupling c is nicely factorized. The Cn\l) satisfy the following W^o algebra 

[£W(l),^likl)] = (n-m)£W^(l) (2-27a) 

[421(1),£W(1)] = (n - 2m)SXma) + mim + 1)/:™ „(1) (2.27b) 

[/:P](i),/:M(l)] = 2(n - m)£|fi„(l) - (n - m)(n + m + 3)£L'U(1) (2.27c) 

and in general 

[£H(1), £W(1)] = i^n - rm)C^:Z-'\l) + ..., (2.28) 
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for r, s > 1; n > —r,m > —s. Here dots denote lower than r + s — 1 rank operators. We 
notice that this Woo algebra is not simple, as it contains a Virasoro subalgebra spanned by the 
£k^'(l)'s. For this reason, it is often called M^i+oo algebra. We also see that once we know these 
generators and £_2(l)i remaining ones are produced by the algebra itself. 

The algebra of the £n' (2) is just a copy of the above one, and the algebra satisfied by the 

[rl ~ [rl 

Wn and by Wn is isomorphic to both. 

The derivation of the W constraints is very simple [^. It consists of taking the coupling 
conditions ( |2.13| ), multiplying them by powers of (5(1) and (5(2), taking the finite trace and 



using the flow equations of the Toda lattice hierarchy. This was done in detail in ref.[15|. There 



one can also find explicit expressions of the generators, see also |22 



2.3 Homogeneity and genus expansion. 

The CF's we compute are genus expanded. The genus expansion is strictly connected with the 
homogeneity properties of the CF's. As we will see the contribution pertinent to any genus is 
a homogeneous function of the couplings (and N) with respect to appropriate degrees assigned 
to all the involved quantities. Precisely, we assign to the couplings the following degrees 

deg{) = [], [ta,k] = y-yak, [N] = y, [c]=y-yi-y2 (2.29) 

where y,yi,y2 are arbitrary constants. Here and in the following N is treated as a coupling 

*i,o = ^2,0- 

If we rescale the couplings as follows 

ta,k A^ "■'^^^Ojfc 

on the basis of the analysis of Bessis, Itzykson and Zuber, |17], we expect the free energy to 
scale like 

oo 

F^Y. A^^^^i-'^^fW (2.30) 

h=0 

where F^^^ is the genus h contribution. In other words 

[fW] = 2y(l - h) (2.31) 

The CF's will be expanded accordingly. Such expectation, based on a path integral analysis, 
remains true in our setup due to the fact that the homogeneity properties carry over to the 
Toda lattice hierarchy. To this end we have simply to consider a genus expansion for all the 
coordinate fields that appear in Q(l) and (5(2), see (|]2l|, |^). The Toda lattice hierarchy 



splits accordingly. In genus the following assignments 

[a^] = (1 + [b?^] = + 1)^2, =yi + y2 (2.32) 

correspond exactly to the assig nments (|]29|) and [F^^)] = 2y. 

It is very common to replace the matrix size with a continuum variable, say x. This is 



permitted provided one rescales all the quantities involved according to the above degrees, Pq] . 
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3 Correlation functions in (unconstrained) matrix mod- 
els. 

We have (at least) three methods to calculate CF's. The first consists of directly solving the W 
constraints; the second consists of determining from the coupling conditions the explicit form of 
(5(1) and Q{2) and then using the flows of the discrete Toda lattice hierarchy; the third method 
is based on passing from the discrete hierarchy to a purely differential one and integrating the 



flows of the latter. The first method has been shown in a number of examples, [l^ and |22]. 
Moreover we will see it at work in the constrained models. Therefore we skip it here and pass 
directly to the second method. 



3.1 Solving the coupling conditions: A42 2 model 



This method is based on explicitly solving the coupling constraints ( p. 13 ). It is then elemen- 



tary to compute correlators by means of eq.( |2.17D and the lattice Toda flow equations. First we 
discuss in detail the model A^2,2; i-e. the purely Gaussian case. For simplicity we set, in the 
following, 

tl,k = ik, *2,fe = Sk 

Lemma 3.1 The matrices Q{1) and Q{2) relevant for the model 7^2,2 fl^^e specified, in 
reference to the coordinatization ( 2.2\ , '2^2^ , by the following coordinates 



/X csi - 2s2ti , . 2s2n 
ao(n) = ao = — — ^, ai(n) = -— — ^, (3.1) 



cti — 2t2Si 2t2n nc 



bo{n) = bo = - — 5", bi{n) = -- — ^, R{n 



452*2 - ' 4s2t2 - ' 4^2*2 " 

The remaining coordinates vanish. 

Proof. The coupling conditions ( p.l3[) for the model 7^2,2 are 

P°(l) + ti + 2t2Q(l) + cQ(2) = 
V°{2) + S1 + 2s2Q{2) + cQ(l) = 

Using the fact that -P^„_i(i) = n for i = 1,2, they can be explicitly written in terms of the 
coordinates as follows. 

2t2R{n) + cbi{n) = cai(n) + 2s2R{n) = 
ti + 2t2ao{n) + cbo{n) = si + cao{n) + 2s26o(^) = 
n + 2t2ai{n) + cR{n) = cR{n) + 252^1 (n) + n = 

These equations can be easily solved and give (^). 

Proposition 3.2 The exact one-point and two-point correlators of the model M.2^2 o,re given 
by the following formulas 



r I 



21=0 k=0 



r!2-^ I N \ / 2^2 \i-k^ 2s2ti-csi Y-2i 

(r - 2l)lkl{l -ky.\ l-k+1 ] Vc2 - 452*2^ V c2 - 4^2*2 ^ 
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Proof. To start with it is convenient to rewrite Q{1) and Q{2) as 

Q{i) = ail+ + (3il + -fie-, i = 1,2 

where 



(3.3) 



+ — ^ En,n+1, ~ X/ ^n,n, ^- ~ X/ ""^".n.-l 

n=0 i=0 n=l 



and 

ai = 1, /3i = ao, 71 
By means of the formulas 



2S2 2*2 . , c 
— — -, a2 = , 132 = bo, 12 = —. 2 



n! 



n=k 



{n-ky. 



N 
k + 1 



we can now make exphcit computations. For example 



Tr(Q(ir) = 



21=0 
r I 

EE 



r 
21 



Tr(/+ + 7ie-)''/3r'' 



-12-^ 



N 



21-- 



^^f^^{r-2l)\k\{l-k)\ \ l-k 



+ 1 



(3.4) 



(3.5) 



(3.6) 



Prom this formula, using ( 2.17| ), we can immediately get eq.( p.2[ ) above. In a similar way we 
can derive < ar >■ 

Pinally using the genus expansion (2.30) we can extract the genus by genus correlators. 
Corollary 3.1 The genus h contributions to the one- and two-point CP's in the model Ai 2, 2 



are 



^ ' {-lf^-'^2~^r\P2h-kil - fc)7V'+i-2^ 

< Tr >h= 2^ 2^ 



21=0 k=0 



(r - 2l)lkl{l - k)\{l -k + l)\ Vc2 - 452*2^ ^ c2 - 452*2 



2S2 \l f2s2tl- CSi\r-2l 



;)( 



(3.7) 



where 



l^k{r) = E ^1^2 ■■■rk, l< k<r, /?o(r) = 1, (3k{r) = otherwise 

l<ri<r2...<rii<r 



Due to the P factor the sum over / in ( |3.7D starts at 2h and the sum over k ends at 2h. 
For the two-point correlators, see Appendix. 

We have given the above proof in some detail since it constitutes a model for all the other 
more complicated cases. In fact there is nothing new when we consider the A^p.i models. They 
can be solved exactly in the same way. New features appear in the case of the A^pi,p2 niodels 
with pi,p2 > 1 and pi + P2 > 4, since the constraints give rise to non-linear algebraic relations 
for the coordinates. Let us see the simplest possible example of this situation. 
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3.2 Solving the coupling conditions: A^3,2 model 

The coupling conditions of the A^3,2 niodel are 

P°(.l) + ^hQ{lf + 2t2Q{l) +ti + cQ{l) = 

+ 2s2Q{2) + S1 + cQ{l) = (3.8) 

Using the coordinatization ( 2.21 ) and ( |2.22| ) we find that the fields ai{n), bi{n), R{n) must satisfy 
the equations 

C62H + ^hR{n)R{n - 1) = 

2t2R{n) + chi{n) + 3t3R{n)(ao{n) + ao(n - 1)) = 
3ts[ao{n)'^ + ai(n) + ai(n + 1) ) + 2t2ao{n) + ti + c6o(n) = 



n + 3t3ai(n) (^ao(n) + ao(n — 1) j + 2^2al(^^) + cR{n) = 

2s2i?(n) + cai(n) = (3.9) 

2s26o("-) + si + cao(n) = 

n + 2s26i(n) + cR{n) = 

One easily realizes that the second, fourth, fifth and seventh equations are linearly dependent. 
Finally one has 

I \ ^ , . ^ Si +cao(n) 

ai(n) = R{n), 60 (n) = 

c 2S2 

, / X n + cR{n) , , , 3*3 , 

^1 " = ^2 n = ^i? n i? n - 1 

2S2 c 



and the recursion relations 

ao(n) + ao(n-l) = -g + -^(c+-^) (3.10) 



2*2 ( 

3t^ ^ 652^3 ^ ^) 

i?(n + 1) + R{n) = — ao(n)2 + f ^ - ^^)ao(n) - + ^ (3.11) 

^ ' ^ ' 2S2 ^ ' V6s2t3 12sit3^ 12^1*3 652*3 

These recursion relations can be solved exactly in complete generality, although the final for- 
mulas may look very cumbersome. However, for our present purposes, it will be sufficient to 
see the solutions in genus 0. In genus the above equations become: 

t \ J,. X , . X si + cao{n) 

ai{n) = R{n), bo{n) - 



c ^ " ' 2S2 

n + cR{n) 3*3 2 /o 

6i(n) = , b2{n) = R{ny (3.12) 

2S2 c 



and the recursion relations 

2*2 c / n 



= + _ _|_)<,„(„, - ^ + ^ (3.4, 



This leads to a cubic equations for ap. Once this equation is solved with the standard formulae, 
all the fields are completely determined. Since they are not particularly illuminating, we do 
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not write down here the explicit solutions. We notice however that, once we know them, it is 
possible to write down immediately an integral expression for the correlation functions. For 
example, using the same formulas as in the previous subsection, one gets 

l>r/2 I I \ I 

where ao, ai are the solutions of the above algebraic system, and we have promoted n and N to 
continuum variables and called the latter x. In a similar way we can obtain all the correlators 
we wish. 

As we see from this example, the method is the same as in the A^2,2 model, the only 
additional difficulty being the solution of a third order algebraic equation. When we come 
to more complicated A^pi,p2 models, the method remain the same but we are faced with the 
problem of solving higher order algebraic equations. 

It has been shown in [^] (see also ||2^) that the model A^o,o represents the c = 1 string 
theory at the self-dual point. Any A^pi,p2 model represents the perturbation of the former by 
the corresponding tachyonic states. We have shown that these perturbations can, at least in 
principle, be solved. However we do not see any point in pushing the analysis further in this 
direction. We would rather like to concentrate from now on on a related interesting problem: 
can one obtain from non-Gaussian Mpi,p2 models 'simple' submodels, in the sense that, for 
example, the correlators are polynomials of the couplings? The answer is yes, and the submodels 
are obtained by imposing constraints in the coordinates of the A^pi,p2 models. The submodels 
are called constrained two-matrix models and to the analysis of some of them are devoted the 
next three sections. 



4 Differential Hierarchies of Two— Matrix Models. 

One possible characterization of the constrained models is by means of the differential integrable 
hierarchies they correspond to. 

Let us return to section 2. We saw there that two-matrix models can be represented by 
means of coupled discrete linear systems, whose consistency conditions give rise to the Toda 



lattice hierarchy. Here we review the method, used in |15], to transform the discrete linear 
systems into equivalent differential systems whose consistency conditions are purely differential 
hierarchies. This is tantamount to separating the N dependence from the dependence on the 
couplings. This section is introduced for the sake of completeness: we collect and try to render 
as plausible as possible the results obtained in |jl5[,|^^ we will need in the following. 

The clue to the construction are the first flows, i.e. the ti^i and t2,i flows. For the sake of 
simplicity let us consider the system I and the flow (|2.24a ) . Let us consider the generic situation 



in which Q{1) has m-i = P2 — 1 lower diagonal lines (see the parametrization ( 2.21[) ). To begin 
with we notice that 

d 



dti. 



= + ao(n)^'n (4.1) 



and adopt for any function f{t) the convention /' = = df. We can rewrite 

= Bn^n+l (4.2) 

where 

-, oo 
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In so doing we implicitly understand that the framework in which we operate is that of the 
pseudodifferential calculus, see for example [p^]. 

It is now an easy exercise to prove that the discrete spectral equation 

is transformed into the pseudodifferential one 

K{l)^n = Al^-n (4.4) 

where 

mi 



L„(l) = d + J2 ai{n)B^_iBr,_i+i . . . B^-i (4.5) 
1=1 

mi 

= 9 + ^a;(r 



1=1 

1 1 



1=1 



d — ao{n — I) d — ao{n — I + 1) d — ao{n — 1) 



Proceeding in the same way for the other equations of system I we obtain the new system 
in differential form 

( \ (^-6) 

4;:*. = (LUl))^^n, 

The subscript + appended to a pseudo-differential operator represents the purely differential 
part of it. The subscript - represents the complementary part. 

Let us come now to the n dependence of the above equations. The operator Ln(l) in 
( [4. 51 ) depends on the coordinates of different lattice points. To deal with this complication, we 
introduce m\ "fields" 5i, ...,5^1, related to the "field" oq in the following way 

Si{n) = ao{n - i) (4.7) 

Then we can rewrite Ln{l) in the following way 

^ d-Si{n) 9-5/_i(n) d - Si{n) 

with the result that Ln(l) is expressed in terms of fields evaluated at the same lattice point. Of 
course the fields St are not independent. However we will consider these fields as completely in- 
dependent from one another in all the intermediate steps of our calculations and only eventually 
impose the condition ([4.7|) . 

To further simplify the notation we will consider henceforth the lattice label n on the same 
footing as the couplings and write 

Oi(n, ...) = ai{n){...), Si{n, ...) = 5j(n)(...) 

where the dots denote the dependence on ti^fci ^2,fc and g. So the expression of Li{n) gets further 
simplified to 

""11 1 
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where, for simplicity, we have dropped the label (1) too. This simplified form is the one we 
constantly refer to throughout this and the following section. 

A similar treatment can be applied to the second linear system as well. Therefore the 
information concerning matrix models can be stored in two differential linear systems + the 
first flow equations ( 2.24c , 2.24%2.24-C , 2.24a )■ The former determine the dependence on the 



couplings, while the latter fix the dependence on N. Therefore what we have accomplished so 
far is the separation of the dependence on N from the dependence on the couplings. 
From now on we refer to the consistency conditions 

|-L=[(L'-)+,L] (4.10) 

where tr = ti,r- ( [1.10| ) are integrable hierarchies, [26|[27|, which are classified by the number 
2m of fields. The pseudo differential operator L in (O) is the relevant Lax operator. 

We can easily locate these hierarchies in a well- 
a particular realization of the KP operator 



cnown framework. In fact L is nothing but 



Lkp = d + Yl 



1=1 



In general, wi are unrestricted coordinates, while in the realizations (O) they are precise 
functions of the fields ai and 5/ and their derivatives. But that is not all, for one can obtain 
new integrable hierarchies via hamiltonian reduction. Each integrable hierarchy characterizes a 
different model. In the case of a reduced hierarchy, we call the corresponding model a reduced 
model. These reduced models will be shown to essentially coincide with the constrained ones. 
The results can be synthesized as follows. 

Summary. Starting from the Lax operator ( [4. 91 ) with given m one can show that: 

1) there are m + 1 distinct differential integrable hierarchies which are obtained by suppressing 
successively the fields Si with the Dirac procedure; 

2) of each such hierarchy it is possible to write down the relevant Lax operator; 

3) at the end of this cascade procedure we find the (m + l)-th KdV hierarchy. 

Therefore for every p = m + \ we have p systems or hierarchies, denoted henceforth with 
the symbol S^, where I counts the number of nonvanishing S fields, < I < m. In particular 
the case I = m corresponds to the 2m-field representation of the KP hierarchy, while I = 
corresponds to the p-KdV hierarchy. 

The above is general and holds for the more complex systems with m > 2. The general case 



was treated in |Q (see also [24]) 



4.1 Examples: the KdV and Boussinesq hierarchies 

The simplest example of L, ( [4. 91 ), corresponds to m = 1. It gives rise to the NLS hierarchy. 

L = d + ai-^ (4.11) 
o — Jl 

If we impose the constraint 5*1 = 0, the second Poisson structure can be reduced via the Dirac 
procedure and leads to the a classical version of the Virasoro algebra. The corresponding 
integrable hierarchy is the KdV hierarchy. Later on we will need the recursion relations for the 
fiows of this hierarchy. They are introduced as follows. Let 

i^2.(x) = 1^ (4.12) 
da[x) 
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where H2k are the Hamiltonians, whose exphcit form can be derived from the Lax operator, 
[^],[^], and a = oi. Then, imposing the compatibihty between the two Poisson brackets, 
characteristic of the hierarchy, we find the recursion relation for the flows 

F2k+2 = DFF2k, DF = d^ + 4ad + 2a' (4.13) 



dt2k+i 
with F2 = a. 

The simplest integrable system that appears in matrix models after the NLS system is the 
four-field representation of the KP hierarchy (m = 2). It naturally leads, via reduction, to the 
Boussinesq hierarchy. Let us describe the latter as concisely as possible. The Boussinesq system 
is described by two fields oi and 02. The Lax operator is 

LB = d^ + aid + a2 (4.14) 

The second Poisson structure is nothing but the classical W3 algebra. The second flow equations 
are 

d d ti 2 

—ai = 24 - a'l, —02 = 4' - -(aia'i + a'/') (4.15) 

This is known as the Boussinesq equation (in parametric form) and it is the first of an integrable 
hierarchy of equations (the Boussinesq or 3~KdV or ^3 hierarchy). 

Like in the KdV case, we give the recursion relations that allow us to calculate all the flows. 
Let us define 

Fr{x) = ^^,\ , Gr{x) = ^^^,\ , ri^Zn 
oa\\x) oa2[x) 

Then imposing the compatibility of the two relevant Poisson brackets, we find the following 
recursion relations 



dai 
dtr 



^=3G;+3 = DcGGr+DGpFr (4.16) 



da2 



3F' = DpcGr + DpFFr (4.17) 



with Fi = 1, Gi = and F2 = 0, G2 = 1. The differential operators are 

Dgg = 3a29 + 24 - aid'^ - 2a[d - < - d'^ (4.18a) 

Dgf = + 2aid + a[ (4.18b) 

Dfg = ^a'^d + 02 - - (aia[ + ajd + 2aid^ + a'" + 3a'(d + Sa'^d"^ + d^) (4.18c) 

DFF = d^ + aid^ + a2 + 3a2d (4.18d) 



5 Correlators in reduced models 

In this section we show that, starting from the p-KdV or Sp hierarchy and borrowing some 
information from matrix models, we can define models, i.e. we can define (and compute) a full 
set of correlators - which turn out to essentially coincide with the constrained models we will 
meet later on. Since however the construction in this section is somewhat heuristic, and, in 
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particular, it does not permit us to carefully fix all the normalizations, we prefer to distinguish 
these models from the constrained two-matrix models of the following section: as they are 
based on reduced hierarchies, we refer to them as reduced models. We call Mp the reduced 
model based on the p-KdV hierarchy. In presenting the reduced models before the constrained 
two-matrix models, which are the true objectives of our research, we follow a historical and, 
hopefully, pedagogical order. 

The essential definition of the reduced models is as follows: we define the correlators by 
identifying the field ai with the two-point function < tiTi >, i.e. we borrow from the matrix 
models eq.( 2.23|) , and differentiate (or integrate) ai with respect to the couplings, as necessary. 



Moreover we only consider the dependence on the tj. = ti,fc couplings and disregard the remaining 
ones. At this point the flows of the relevant hierarchy allow us to calculate the correlators, at 
least up to some constants - we are going to see some examples later on. The reason for this 
is as follows. A part of the information concerning the coupling conditions is in fact stored in 
the differential system of the model: the Lax operator inherits the information contained in 



the second equation (2.13) via the number of non-vanishing diagonal lines of the original Q{1) 
matrix. Therefore it is not surprising that the flow equations are almost enough to determine 
the CF's. However not all the information concerning the coupling conditions is contained in 



the differential hierarchy which characterizes the model, the first equation ( 2.13| ) is not, and this 



is reflected in the undetermined constants that appear when we try to calculate the correlators. 
Let us see this point in detail in an explicit example. 

5.1 The KdV hierarchy and the associated A^2 model 

We showed in section 4 that we are allowed to impose the constraint S = on the NLS system 
while preserving integrability. In other words there is a consistent subsystem of the NLS system, 
of which we can easily compute the flows, ( 4.13| ). These are the KdV flows. We recall that only 



the odd flows survive the reduction. Therefore the t2n are disregarded. It is therefore natural 
to forget to = N as well. 

To start with let us define the critical points for this model: 

The k-th critical point of the M.2 model are defined by 2{2k + l)t2fe+i = —1 o,nd ti = for 
lj^l,2k + l. 

For the origin of this terminology and further properties of critical points in matrix models 
see |p2|| . We will see next that the above critical point corresponds to a two-matrix model with 
a Vi potential of order 2fc + 1 and a V2 potential of cubic order. Let us notice 1) that the 
correlators of Ai2 at the various critical points are functions of ti alone, 2) that in order to 
preserve the homogeneity properties at the /c-th critical point we have to set y = yi{2k + 1) in 

In the following we study the first critical point, k = 1. On the basis of the assignments 
of section 2.3 we have [a] = [ti] = 2yi, in genus 0. Therefore it must be: a ~ ti. The 
proportionality constant can be absorbed with a rescaling (provided it is non-vanishing, which 
is the case as we shall see). So we start from the position 

a = ti (5.19) 

Then we have: 

Lemma 5.1 As a consequence of ( 5.H ), the functions i*2n relevant to the KdV hierarchy 
are given by 

00 

F2n = a{n, h)t"'-^^, n>3h (5.20) 

h=0 
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where 

a{n, h) = ^— 7 V- 7TT, (5.21) 

and rW is the 2 by 2 factorial, i.e. r\\ = r{r — 2) . . . up to 1 oi 2. 

Proof. We insert the expression (|5.20|) into (|4.13D and obtain the recursion relation 



2n — 6/i + 1 

a(n +l,h) = {n-3h + 2){n - 3/i + 3)a(n, h - 1) + 2 a(n, /i) (5.22) 

~|~ 1 



for the coefficients a{n,h). One can immediately verify that (5.21) is a solution of ( |5.22| ), but 



it is not unique. While integrating ( ^.22 ) one has to specify what = a{3h,h) are V/i. The 



latter are the coefficients of [a')"^^ in Fg/i ^ind satisfy the recurrence relation 
3hbh = 2{6h - l)(6/i - 3)(6/i - 5)bh-i, 



One immediately gets ( 5.21[ ) . This ends the proof of the Lemma 



Proposition 5.1 The exact one-point correlators of the A^2 model at the first critical point 



are 

in— 1 



where the genus expansion is explicitly exhibited. 

Proof (partial). We have simply to recall that =< T2n-iTi > and integrate over ti. We 
obtain ( |5.23| ) for n > 3h. The values of < rg/i-s >h (i-e. n = 3h — 1), which are obtained by 
simple extension of this result, are also correct, but strictly speaking they do not follow from 
the previous argument: they are pure integration constants and cannot be obtained from the 
flows alone. We will be able to completely justify eq.( 5.23| ) only by means of the W constraints. 



It is in fact the VF-constraints that completely determine such constants. 

5.2 VF— constraints of the A^2 model 

Some information concerning the coupling conditions is not contained in the differential KdV 
hierarchy. In order to retrieve it we have to use the W constraints. The point is that they cannot 
be the same as in the unreduced models, since the hierarchy underlying the model has changed, 
and we recall once again that the W constraints are based on the flow equations. Therefore we 
have to reconstruct effective W constraints on the basis of the reduced hierarchy. Let us argue 
as follows. In reduced models we are interested in solutions that do not depend on the second 
sector (i.e. on t2,fc)- If we look at eq.(|2^), we see that such solutions should therefore satisfy 



to VF-constraints of the form 

£M(l)ZAr = 0, r>l; n > -r (5.24) 

Consequently we look for W constraints of this type, with generators belonging to a algebra, 
which are however compatible with the KdV hierarchy. It is easy to see that the universal 
generators in ( 2.25 ), (see [^]), do not satisfy the KdV flows. We find instead 



Proposition 5.2 The effective W constraints for Ai2 take the form 

LnVz = 0, n>-l (5.25) 
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where 



^0 = + 9)^7^7 + 777 (5.26) 

in = E(^+ 9)^2^+1 7^7 +E7^7 JTf ' 

^ 0't2fe+2n+l j^,^Q Ot2k+lOt2n-2k-l 

These generators satisfy the commutation relations of the Virasoro algebra. 

Proof. Let us prove first that ( ^.25 ) are in agreement with the KdV flows. To this end we 



differentiate (5.25) with n > with respect to ti. Using the definition of F2k, we can write 



(remember the notations introduced after eq.(4.1)) 



^(2A: + l)t2fc+ l-^2fc+2n+2 + 9 F2n+2 



k=0 

n-1 



+ (277;^ F2n-2k + F2k+2d ^ F2n-2k + 9 ^-F2fc+2-^2n-2fc 



Here is understood in the sense of the pseudodifferential calculus and denotes indefinite 
integration (see below for further specifications). Now we apply to it the recursion operator Dp. 
What we obtain, by using eq.( [4.13| ), is nothing but the constraint Ln+i^/Z = differentiated 
twice with respect to x. To see this we have to apply the remarkable formula 



n-l 

F2n+4 = F2n+2 + 3^2-^21+2 + ^ (2^2^+2-^21-2/0 " ^2A:+2^2n-2fc 

fc=0 

+4F2F2fc+2-^2n-2fc — -^2fc+2-f2n-2fc+2 



which can be obtained once again from the recursion relation ( 4.13 ). As for the cases n = and 
n = —1, which have not been included in the above argument, they can be explicitly verified. 

What we have done so far amounts to saying that starting from L-\\fZ = and successively 
applying the operator O = d~^Dpd, we obtain all the Ln\fZ = 0. Here we have to exercise some 
care with the double integration d^^. d"^ represents an indefinite integration which preserves 
the homogeneity properties. This is a perfectly well defined operation unless the output of it has 
degree 0. In such numerical integration constant may appear. Now, in Z-i/2l„z1/2 

there appear contributions of degree 2y(l — h) + 2nyi, with n = 0, 1, 2, . . .. So, since y and yi 
are generic numbers, the only dangerous case (in the above sense) is when h = l,n = 0. In 
other words when we pass from L_i\/Z = to Lq^/Z = by applying O we are not guaranteed 
that the appropriate constant is given by the jq present in Lq. However at this point we make 
the request that = 2Lo, and this unambiguously fixes such constant. It remains for us 

to justify -L_i\/Z = (which is often referred to as the string equation) or rather the term ~ tf 
in From the degree analysis one sees that the only possible polynomial of the couplings 

one can write is t^. Its coefficient is determined by the requirement that, applying the recursion 
device to one gets Lq. 

Finally, we do not look for higher tensor constraints since eq.( 5.25| ) is enough to determine 
everything. 

On the basis of eq.( p.25| ) one can complete the proof of Proposition 5.1. 
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5.3 The Boussinesq hierarchy and the associated A4s model. 

The 3-KdV or Boussinesq hierarchy was introduced in 4.2 as a reduced hierarchy. The corre- 
sponding model is denoted 7W3. It is described by two fields ai and 02 and is specified by the 
Lax operator 

L = d^ + aid + a2 (5.27) 
In the Boussinesq hierarchy the t^k flows with A; = 0, 1, 2, 3... do not appear. 



The correlation function interpretation of the fields ai and 02 is given by eq.( 2.23 ) and the 
first of ( p5D : 

ai =< TiTi >, 2a2 =< T1T2 > + < TiTiTi >, (5.28) 

Now we proceed as in the KdV case (but skip many details). The first critical point is 
defined by 

4*4 = -1, tk = k> 4 

This implies in particular that y = 4yi and that the correlators will be functions of ti, ^2- Next, 
the degree analysis shows that ai ~ t2 and 02 ~ ti. An elementary use of the first eq.(p^ 
shows that, up to an overall multiplicative normalization constant, we can choose 

ai = 6*2, a2 = 3ti (5.29) 

This will be our choice (as it is consistent with the W constraints and the definition of the 
critical point). Now it is relatively easy to use the recursion relations of the flow equations to 
compute the correlation functions. 

Proposition 5.3 The exact one-point correlators of M.^ are 

00 

< T"3n-2+e > = X/ ^ >h 
h=0 

< r3n-2+. >h - 1^ 48^^/1! 3"-i-2+E(^_i + g)! (3n - 3i + 2 + e - 8/1)!! ' 

n-i+ee2Z 

(3n-3 + 2e)!!! (3n-2 + 2e)!!! 3n-3.,+2+. _^^ 



(3i)!!! ( 3n~3.R. _3^^,,, ( 3n-3.j+2+. _3^^,,, 



where e = 0, 1 and n\\\ is the 3 by 3 factorial, i.e. n!!! = n(n — 3)(n — 6) . . . up to either 
1 or 2 or 3. By convention 0!!! = (—1)!!! = 1, and ^ ~ ~ ^ f'^^ ^ — ^ consequence 

in the above formula the exponents of ti and t2 are always non-negative. 

Sketch of proof. One has to remark first that, while the contributions from two contiguous 
genera differ by 8yi, the recursion operators Dqg,---., Dff contain contributions that differ by 
4yi. It follows that in Gr, there will appear half-genus contributions. Therefore, at the first 
critical point, we have to start from the ansatz (case e = 0) 

3n — 3j — 1 



G3„+i = G^zll\= E o.,{n,s/2)t{t7'^-''^-^' (5.31) 

n-je2Z+l 

W = E^3nS> Fti\= E fiM^s/2)t{tf^-'' (5.32) 



3n — 3j 



-2s 



s=0 j=0 

n-je2Z 
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where s is the half-genus label, i.e. s = 2h, and the exponents of ti and t2 are always non- 
negative. The half-genus contributions must not appear in the correlators, thus we must have 
the physicality conditions 



aj{n, s/2) = 0, 2/3, (n, s/2) = (j + l)a,+i(n, 



s G 2Z + 1 



(5.33) 



Plugging ( ^.32 ) and (5.31) into (4.16,4.17), and using ( |5.33| ) we find the following recursion 
relations for the coefficients 



3jaj(n + 1, /i) = 3(3j - l)a,_i(n, h) + 12j/3, (n, h) + 2j{j + + 2)/3,+2(n, h - I) (5.34) 
for n — j £ 2Z, and 



3j(3j{n + l,h) = 3(3j - 2)/3j_i(n, h) - 2Aja,{n, h) - 5j{j + + 2)aj+2{n, h - 1) 



j{j + l)(j + 2)(j + 3)(j + 4)a,+4(n, h - 2) 



(5.35) 



for n — J G 2Z + 1. These relations can be integrated and give the following result: 



aj{n, h) 



Zn — Zj — l 

1 2 — 2 



■3/i 



(-1) 



(3n-3j-l-6/i)! 



48^/i! 3"-i-i(n-l)! (3n - 3i - 1 - 8/i)!! 
(3n-3)!!! (3n-2)!!! 



(3j)!!! ( ^"~2^''^ - 3/i)!!! (3!!^3i^ - 3/1)!!! ^ 



n - j e 2Z + 1 



and 



Pj{n,h) 



3n — 3j 

1 2^' 



3h 



(3n-3j-6/i)!! 



AS'^hl 3"-J-i(n-l)! (3n - 3i - 8/i)!! 

(3n-3)!!! (3n-2)!!! 

' (3j)!!! (2^ - 3/i)!!! (^^ - 3/i - 2)!!! ' 



n-je2Z 



Now we recall that 



< Trn >h= 3G|,'23, < rrT2 >h= 6F, 



r+3 



and integrate w.r.t. ti and t2 the first and second expressions, respectively. Comparing the 
results we find (153C| ). Just as in the KdV case we have to treat separately the case when both 
exponents of ti and ^2 in (5.30) vanish. The coefficients given by ( p. 301 ) are the correct ones, 
but they have to be checked by means of the W constraints. 
Likewise we can compute < r^n-i > (case e = 1^. 

The effective W constraints in the case of the model are found once again by requiring 
that they be consistent with the Boussinesq flows and that the W generators form a closed 
algebra. 

Proposition 5.4. The effective W constraints for the ^A■i model are 



W Z 



1/3 



0, 



1,2, n>-r 



(5.36) 



where 



00 

-Ykt 



fc+i=— 3n 



Vn 



k.l 
k+l=3n 
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^l,,l,=l Mh+l2+3n 9 fj^. dtkdtj 

l—k—j=—3n 

1 9^ 1 

i+fe+j=3n ;+fe+j=-3n 

in i/iese expressions summations are limited to the terms such that no index involved is either 
negative or multiple of 3. The above two sets of generators form a closed algebra, the W3 algebra, 



n+m 



[lL^UI^I] = -^(n-m)((n2 + m2 + 4nm) + 3(n + m) + 2)4i„ 



where 



k<-l k>0 

This corresponds to the quantum W3 algebra with central charge 2. 

Sketch of proof. One can prove the consistency of the Boussinesq flows with the above 
constraints in the following way. Call Kn = Z~^/^Ln^Z^/^. Then, for example, one can check 
that 

and so on. In fact it is not necessary to prove such equalities for any n > —1, we simply 
need to do it for the first few cases, for Kn = for n = —1,0, 1,2 implies = 0,Vn > — 1 
due to the Virasoro algebra structure. The same can be done for the higher order constraints. 
Constants and polynomials in the couplings which appear in the generators can be fixed by 
simply requiring algebraic closure. 



5.4 Other models 

Let us generalize what we have just done for the Boussinesq hierarchy to the Aip (or p-th 
KdV) models. The general recipe is as follows. One must first of all disregard all the tk with k 
a multiple of p; the first critical point is 

{p + l)tp+i = b, tk = k>p + l (5.37) 

where b is any number. The degree assignment is 

[tk] =p+l-k, [fW] = {2p + 2)(1 - h), [aS°)] = 2, . . . , [afl,] = p (5.38) 

where af^ is the genus part of Oj and we have set, for simplicity yi = 1. The CF's will be 
homogeneous functions of ti, . . . , tp-i, which constitute the small phase space. 

In all the cases the method to compute CF's is the same as before. We do not have however 
to redo literally the same steps as before. A shortcut consists of fixing the form of the fields 
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by means of effective W^-constraints, which in turn are determined imposing compatibility with 
the relevant flow equations. Once this is done the CF's can be obtained from the flow equations. 
We write down hereafter the -constraints for the general A4p model. The W constraints 

are 

LMz^/f = 0, r=l,...p-l, n>-r (5.39) 

Compact formulas for the above generators can be written down by means of the bosonic 
formalism. Let us introduce the current 

OO -, OO o 

r=l ^ r=l 

Then 

pT-^ \ J r + 1 

The normal ordering in the last definition is the one between ty and The derivative is 
always supposed to stay at the right. These generators close over a Wp algebra with central 
charge p — 1. In particular we have 

In the above formulas n is any integer, and multiples of p as well as non-positive integers are 
excluded among the summation indices. 

We can extract particular exact formulas as follows: we write down the dispersionless version 
of the constraints L^^j.Z^/p = 0, with r = 1, . . . ,p — 1; this equation gives a recursion relation 
for < >, r = 1, . . . ,p — 1, in terms of < > with / < r, which can be solved and gives: 

J — 1 

p-l n-l / \ 

<rr> = E (-1)" n E b-- (2) E 4^^ • • • • • • t,,, ■ 

n=l i=l Si=l \ * / (i) (i) ' 

'i '■■■•'■Si 

l['^+...+li[Un-s,(p+l)=ri+i 

■ E b-^^'- {"t^] E 4"^ • • • • • V") (5-40) 

n \ " / , ; , ^ 1 r-s„ + l 

Sn=0 \ / j{n) j(n) 

1 r — sn + 1 

Although this result has been obtained from the a genus approximation, it is an exact result. 
In particular, setting b = —1, we have 

< n >= I E ''^^kti (5.41) 

k,l;k+l=p 

As a consequence 

< TiTkTi >= pkl5i^p_k, I < k,l < p - 1 (5.42) 
which specifies the metric of the corresponding topological field theory, [p^]. 
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5.5 Higher critical points 

In all the previous examples the first critical point has been characterized by a dependence 
of the basic fields on the couplings specified by homogeneous polynomials with non-negative 
integer powers. Higher critical points are characterized still by a homogeneous dependence, but 
with rational and/or negative powers of the couplings. 

The procedure to compute correlation functions is always the same: a shortcut to arrive 
at the results is to use simultaneously W constraints and flow equations. With a good deal 
of perseverance we could probably arrive at exact correlators as in the previous subsections. 
However, in order to give an idea and for future reference, we think it is enough to present a 
few partial results. 

Let us start with the KdV model. The (Kazakov's) critical points were defined in subsection 
5.1. At these points the degree assignments (setting yi = 1) are 

[ti] = 2k + l-l, [F^'*)] = 2(2A: + 1)(1 - /i), [a^^^] = 2 - {2k + l)2h (5.43) 

Indeed, contrary to the first critical point, we have to expect non-vanishing contributions from 
all genera for the field a = a^^'^ + a^^^ + a^'^^ + .... We find the following results 

a(0) = a{k)'/''ty\ a{k) ^' 



2''-H2k - 1)!! 



4fc+l 

a(2) = -a{k)-^'^-f{k)t^ " 



where 



lik) = ^^(k - l)(32/c3 - 72k^ + Ulk - 77) 

Knowing this formulas one can calculate the correlators in genus 0,1 and 2. The expression for 
a^^^ is also in ref.|2S] 



Higher critical points of more complicated models can be found in p| 



The models A4p have a topological field theory interpretation. The corresponding topological 
field theories are easily identified with those of the A series in the ADE classification The 
latter are known to be based on p-KdV hierarchies. Therefore, what we have achieved in this 
section is a new presentation of this old subject, in fact a very powerful presentation since it has 
allowed us to calculate new all-genus expressions for the correlators. However, if we look at this 
section not from the point of view of topological field theories but from the matrix model point 
of view we cannot yet be satisfied. Although we have used the reduced integrable hierarchies 
obtained in section 4 from two-matrix models, and we have used other matrix model inputs, 
a direct connection between two-matrix models and the results found in this section, although 
very plausible, has not yet been established: in particular some normalization constants have 
been arbitrarily fixed and the identification ai =< tiTi > deserves a safer ground. 

The purpose of the next section is to provide such connection. 

6 Constrained two— matrix models. 

In this section we want to introduce and analyze constrained two-matrix models that are char- 
acterized by p-KdV hierarchies. We are going to study in particular detail the ones based on 
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KdV and Boussinesq hierarchy. In the previous two sections we learned that, in order to end 
up with the p-th KdV hierarchy, we have to impose the constraints Si = 0, which amounts 
in the Toda lattice formalism to impose that the diagonal elements of the matrix Q{1) vanish 
identically. We will just impose this constraint on the model Mpi,p2, and call the two-matrix 
models so obtained A^pi^p2- The main result of this section is that the reduced models Aip can 
be imbedded in the renormalized -A4p+i p- 

A remarkable aspect of our derivation is that we obtain all the results via the flows of the 
Toda lattice hierarchy and we never abandon the framework of two-matrix models. A byproduct 
of this section is a new way to derive p-KdV hierarchies from Toda lattice hierarchies. 

6.1 The constrained two— matrix model ^^3*^2 the KdV hier- 
archy 

The simplest interesting model is M^^l- start from the coupling conditions of the 

model A^3,2 of section 3.2, and impose the condition ao(n) = 0. 

6.1.1 The coupling conditions in 

Remember that KdV, at the first critical point, is expected at t2 = and 6t3 = — 1. Imposing 
this and 

ao(n) = (6.44) 

eqs.(^) becomes, in genus 0, 

2cb2 -R^ = 0, 61 = 0, ai = ti - — 

2S2 

R = , bo{n) = 2s2R + cai = (6.45) 

Setting for simplicity si = 0, one gets in particular 

ai = ti (6.46) 

The importance of ( |6.46| ) is that, from the constrained coupling conditions, we have obtained 
oi = ti, a result which was postulated in section 5.1 on the basis of plausibility arguments. 

Remark 1. The particular values chosen for si,t2)^3 are not important. Other choices 
would not change qualitatively the results, but only renormalize them either additively (si) or 
multiplicatively (^3,^2)- 

Remark 2. The last equation ( |6.45 ) implies 

ti = (6.47) 

The conditions ( |6.46| ) and ( |6.47| ) have to be imposed as a last condition on the correlators after 
all the calculations have been carried out (see below), therefore the fact that ti ~ n does not 
interfere with differentiating or integrating with respect to n. The rather mysterious condition 
( |6.47| ) may have an interesting topological field theory interpretation, see |23]. 



Remark 3. We stress that it is irrelevant whether oi = ti is the unique solution of the 
coupling conditions and that it be found in genus 0. What is important is that we be able to 
impose it (at every genus) without breaking integrability. This is in fact what we are going to 
show next. 



24 



6.2 KdV flows from Toda flows 



We show next that we can extract the KdV flows directly from the Toda lattice hierarchy. 
We recall that in M32 



Q{1) = /+ + ^ [ao{i)Ei^i + ai{i)Eu-i) 



Setting Q = Q{1), we have 

Lemma 6.1. In 7^3,2 the formulas, obtained from the Toda lattice hierarchy and from 
eq.^Hj) 

<rirfc >=Tr([Q+,Q'=]) (6.48) 

give rise to the flows of the NLS hierarchy. 

Proof. We notice that, due to eq.(2.23), we have 

dai d 
Otk oti 

This is the first ingredient. The second ingredient are the first flows ( p. 24a ), which, in the A^3,2 
case, take the form 

S[ = {1- D^^)ai, a'l = ai{D^ - l)Si (6.49) 

(remember that we defined Si{n) = ao(n — 1)). Here we have introduced a new notation, which 
turns out to be very convenient in this kind of trade. For any discrete function /tv we define 

{Dof)N = fN+l 

We will also use the notation e'^° instead of Dq, with the following difference: Dq applies to the 
nearest right neighbour, while is meant to act on whatever is on its right. 
Using this notation we can write Q as follows 

Q = e^o + DqSi + aie-^° (6.50) 

We also remark that the sum, J2n=o^ exactly the inverse of the operation Dq — 1. Now 

it remains for us to evaluate the RHS of ( |6.48D . Hereon we give an example; a more complete 
proof will be provided elsewhere 

dai d 



_d 

dt 



-[Doaiai + oiDq ^ai + aiai + {DoSi)'^ai + Sfai + aiSiDoSi 
1 ^ 

= (^a'( + 3al + SSfai + 3Sia[y 

which is exactly the third NLS flow for ai. In order to flnd the flows for 5i one has simply to 
differentiate 

Si = d'\l - Do')ai 

and use again the first flows. 
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Using this Lemma, it is now easy to conclude our argument. We pass to the Mfl model 
setting = 0. The odd flows of oi reduce to the KdV flows in exactly the same form given by 
the recursion relations ( [4.13| ) . Not only are these relations compatible with the W constraints of 
section 5.2, but also the result oi = ti, obtained from the coupling conditions, is, and coincides 



with our assumption (5.19). We conclude that the model Ai2 is imbedded in the constrained 
two-matrix model and that all the results obtained in section 5.1,5.2 hold true for the 

latter. 

Remark. The second sector of i.e. the dependence on s^, as well as the dependence 

on the bilinear coupling c, does not play a role in the above arguments. They can at most 
renormalize the final results (as noticed above). In other words, the second sector is a spectator. 
Whether and how it is possible to compute correlators of the second sector, i.e. correlators of 
(Tfc, or mixed ones, is a question that we leave open here. 

6.3 The k-th KdV critical point and M^2k+i.2 

To confirm the result just obtained let us look at the k-th KdV critical point (section 5.5). This 
critical point turns out to be embedded in the model A^2A:+i 2- '^^^ relevant potentials in this 
case are 

^l(Al) = WlAf +^ + tlAi, V2{\2) = S2\l + S1A2 

which entail the following coupling conditions 

P°(l) + {2k + l)t2k+iQ{lf^ + ti+ cQ{2) = 0, P°{2) + 2s2Q{2) + si+ cQ{l) = 
The relevant equations one gets in genus are 

{2k + l)t2fc+i ( ] all + c6o + ti = 



252^0 + Si = 0, 2S2R + cai = 0, 252^1 + cR + n = {) 



One finds 




60 = -^, «! = (- 

Here, either we make the replacement ti — > ti = ti — ^ or simply set si = 0. Moreover we set 
(k-th critical point) 2(2A; + l)i2fc+i = -1- Then 

f k\ \i/k 

= i 2^-i(2fc-l)!! ^0 (6.51) 

This result coincides exactly with the analogous formula in section 5.5. Substituting this in the 
flows we can calculate the KdV correlators for any critical point in genus 0. 

6.4 The renormaUzed two— matrix model ^^4*^3 and the Boussi- 
nesq hierarchy 

We want now to repeat the same for the Boussinesq hierarchy. We have to start from the 
but one soon realizes that things are not as simple as in the previous example. In fact an analog 
of Lemma 6.1 holds, but setting oq = in the resulting flows does not lead to the Boussinesq 
flows. A significant change of strategy is necessary. But let us start once again from the coupling 
conditions for ^^4,3. 
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6.4.1 The coupling conditions in A^i^s 

The coupling conditions for the two-matrix model 7^4,3 are 

P°{1) + iUQilf + 3*3^(1)^ + 2t2Q(l) +ti + cQ(l) = 

P2 + 3s3Q(,2f + 2s2Q{2) + si + cQ{l) = (6.52) 

These can be expressed in terms of the fields ao, ai, 02, 601 ^i, &2, ^3 and R. We expect to find the 
Boussinesq hierarchy for ag = and ^3 = 0. The coupling conditions in genus then become 

C63 = -AUR^, cb2 = 0, cbi = -2t2R - 12UaiR, 

cbo = -ti - 12t4a2, n + I2t4al + 2t2ai + cR = 0, Ss^R"^ + ca2 = 0, 

6S360-R + 2s2i? + cai = 0, 53(261 + 6g) + 2S2&0 + si = 0, 65360^1 + 2^261 + ci? + n = 

For simplicity we choose 61 = and si = 0, which leads to 

"^ = "6^' ^^ = -T2tI' ^ = -c ^^-'^^ 

After determining bo = and 63 we are left with two conditions of the type ( 6.47| ) among the 
couplings, which are irrelevant for the following developments. Different choices of 61 and si 
would simply imply additive redefinitions (renormalizations) of ti and t2, therefore we ignore 
them. 

It is important that up to a global rescaling (see below) we have found the same results that 
were assumed under (very) plausible arguments in section 5.3, see eq.( 5.2{\ ). 



6.4.2 Boussinesq flows from Toda flows 

As we have anticipated above, in the Al4,3 there is an analog of Lemma 6.1. In particular, from 
the Toda flows we can obtain the flows of the four-fleld KP hierarchy. However this is irrelevant 
to our problem, since, setting Si = S2 = does not lead to the Boussinesq hierarchy (or to any 
integrable hierarchy, for that matter). The reason is well-known: the above conditions oblige 
the system to flow outside the manifolds of the flow equations. To preserve integrability we have 
to introduce in the original hierarchy a (presumably infinite) set of corrections. In field theory 
language we can say that the constraint = can be imposed without spoiling integrability 
only at the price of introducing a (presumably infinite) set of counterterms. The very important 
point is that this set of counterterms can be exactly computed, after which the resulting model, 
referred to henceforth as the 3 model, will accomodate the Boussinesq hierarchy. 
The recipe to obtain the result is as follows: 

1) Define the general matrix 

00 

Q = + ^ aie-'^° (6.54) 
1=1 

2) Assume as first flows 

Dqcli = ai, Doai = dj + a[_i, i>2 (6.55) 

(these are the usual flrst flows in which we have set do = 0). 

3) Now impose 

q3 ^ gSao _^ 3Q^g9o _^ (6.56) 
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This equation and (6.55) completely determine in terms of ai,a2, 



2 

«i = ai = ai, 02 = ^2 = «2 — a'l, 03 = as = — 03 — + -a'/, 
2 1 

04 = £14 = —02 — —0,1 + 4aia]^ — 2aia2 

and so on. Q is now to replace the matrix Q{1) of the model A^4,3 with ag = 0. It contains the 
right counterterms to generate the Boussinesq hierarchy. To this end, 
4) use the following 

Lemma 6.2 Replace Q{1) with Q in the Toda lattice hierarchy formulas and evaluate them 
at hi = Oi 

|| = ±Tr{[Q+,Q%.=-a.. k + 3n (6.57) 

and the like. Then these formulae provide a realization of the Boussinesq hierarchy. 

Proof. We limit ourselves to a few examples. A more complete proof will be given elsewhere. 

^ = ^'^^(['3+' Q^]) \a^=a, = {Dod2 + a2)'|a,=a. = 2^ - o" 

dai _ d _ 2 „ 1 



= Q ] j = - gfli' - 2aia'i + 4aia2 

da2 1 d „ /r^2 r>2i\\ 1 , 2 2 „ 

and so on. These are the Boussinesq flows, even though not in exactly the same form as in 
section 4.1 and 5.3. In fact we have to multiply by 3 the ai, a2 fields there to obtain the flows 
here. 

Remark 1. We should have allowed also for a field oq in Q and set it to zero at the end of 
the calculations. This could of course have been done but would not have changed the results. 
Setting oq = from the very beginning we have simply anticipated the result and simplified the 
formulas a lot. 

Remark 2. It is not surprising that we have found some disagreements in the normalizations 
here compared sections 4 and 5. We have already pointed out at the end of section 5 that some 
identifications made there were likely to be arbitrary as far as the normalizations are concerned. 

Up to the normalization problem illustrated in the previous remark we can see that the flow 
equations, the coupling conditions l\6. 5S ) and consequently the W constraints pertinent to the 



A^3 model can be embedded in the M\ 3 matrix model. If we set the critical point at 12t4 = — 1 
in ( |6.53| ) and multiply by 3 the fields ai, ^2 in section 4 and 5, we can simply transfer the results 

fO t) 

obtained there to A^4 3 . To be more precise, we summarize the results concerning the latter 
as follows: 

i) The Boussinesq flows are given by 

^ = G;+3 = DGcGr+DGpFr (6.58) 



da 



2 



F;+3 = DpcGr + DpFFr (6.59) 



dtr 

with Fl = 1, Gi = and F2 = 0, G2 = 1. The differential operators are 

Dgg = 3a2d + 2^ - aid'^ - 2a[d - a'{ - \d'' (6.60a) 
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Dgf = i:d^ + 2aid + a[ (6.60b) 

2 1 1 

Dfg = + 4 - 2(aia'i + aid + -aid^ + -a'/' + a'/a + a[d'^ + -9^) (6.60c) 

= la"^ + ai(92 + + 302^ (6.60d) 
3 



a) The coupling conditions imply, at 12t4 = — 1, 

ai = 2*2, a2 = ii (6.61) 

in) The correlators are the same as in section 5.5 except for a global factor of 3. In particular 
the RHS of ( |5.30| ) must be divided by 3. 

iv) The W constraints appropriate for M.^^'^^ are 

lM Z = 0, r = 1, 2, n > -r (6.62) 

the generators being the same as in Proposition 5.4. But in order to reproduce the right 
correlators we have to compute them at 4*4 = — 1. In other words the location of the critical 
point gets renormalized. 



7 Conclusion 

The procedure introduced above for the Boussinesq hierarchy holds for any p-KdV hierarchy. 
We can always define suitable coordinates in the matrix which, inserted in the Toda lattice 
flows, generate the p-KdV flows. The recipe is just a generalization of the one given above. This 
is certainly a remarkable (and new, to our knowledge) result, which deserves further elaboration. 

Finally we can draw the following conclusion: the p-KdV hierarchies are contained in specific 
constrained two-matrix models; once we impose the constraint, the Q matrix of the relevant 
model has to be suitably redefined (except in the 2-KdV case) to insure integr ability; the 
counterterms can be exactly calculated and give rise to "renormalized" coordinates; in turn 
these coordinates, when substituted in the formulas of the Toda lattice hierarchy, gives rise 
to the p-KdV flows. The same procedure may well be applicable to extract from two-matrix 



models other hierarchies such as those studied in [30| 



Appendix 

Here are the exact 2 point CF's of the model A^2,2- Let us first define the functions 



1/2 '-^ m p/2 . . 

n\m\ 



^n,m(a2,^l,^2,7l,72) = Y.Y.Y. ^kT^ 

1=0 k=0 r=0 p=0 0=0 



p=ug=(i ,n - l)\{m - p)\klq\ 



1 



min{l—2k,p—2q) 

fr'o ^r\j\{r + k-q+ -2k- r- i)!(^ - k - q - r - j)\ 

1 

-2k- r)!(r - i)!(^ -k-q- r)!(2^ -q + k + r-j) 

'^+P U „ .\\{ ^ \ J^^^+''- on-lnm-p k^rj^-k- 
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where 02, /32j 71, 72 were defined in eq.(^). Then 

<Tnam> = TY([Q(ir,Q(2n) =.F„,„(a2,/3i,/32,7i,72) (7.63) 
<TnTm> = Tr ( [^(1)", Q(in ) = A , /?! , 7i , 7i ) (7.64) 

Finahy < GnCTm > is obtained from < T„Tm > with the exchange <-> for k = 1,2. 
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